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Going out with a bang: compact object collisions resulting 
from supernovae in binary systems 
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ABSTRACT 

Binary star systems containing a neutron star or a black hole with an evolved, massive 
star are dynamically perturbed when the latter undergoes a supernova explosion. 
It is possible that the natal kick received by the newly-formed neutron star in the 
supernova may place the stellar remnants into a bound, highly eccentric orbit. In 
this case, the two compact objects can tidally interact and spiral into one another 
on a short timescale. The interaction with an accretion disc of supernova debris is 
also considered. We quantify the likelihood of such events and show that they would 
be expected to produce a high-energy transient, possibly a short gamma-ray burst, 
typically within a few days of the supernova. 
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1 INTRODUCTION 



The coalescence of a neutron star - neutron star (NS-NS) 
binary under the action of gravitational wave emission is 
arguably the leading progenitor mod el fo r short gamma-ray 
burst s (GRBs: 
19921: 



e.g. [ Eichler et al.l Il989l: iMeszaros fc Reea 



shocks, and for systems in which one of the neutron stars is 
replaced by a black hole (BH). 

Interest in such collisions has been rekindled by the ob- 
servation of prolonged central engine activity in some short 
GRBs. Though their nature is still under debate, features 



Naravan et all Il992l: iRuffert fc Jankal 1 19971 : 



Rosswog et all l2003l : lOechslin fc Jankal |2006| ). Colli 



such as X-r ay flares |Burrows et al. |2005|) or long soft tails 
of emission (iBarthelmy et al.ll2005l : Worris fc Bonnellll2006l ) 



sions between two neutron stars also do occur, for instance 
in regions of large stell ar number de n sities, but at a 
substantially lower rate. iKatz fc Canell (|l996l ) suggested 
such events might be the central engines of short GRBs. 

In gravitational wave-driven coalescences the extremely 
large sound velocities in nuclear matter, Cs~0.4c, result in 
subsonic initial encounters, and shocks only form in later 
phases after rotationally shed mat erial interacts w i th it- 
self and with the centr al object (|Ruffert fc Jankal I2001I : 
iRosswog fc Davie^ |2002| . in particular their Fig. 13). In 
stark contrast, direct collisions of two neutron stars ap- 
proach each other with relative velocities close to free-fall. 



^'frocfall ~ 0.64 c 
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, and therefore pro- 



^2.8Mq ) \Ri+R2 J 

duce strong shocks pa ssing through the nearly unperturbed 
initial neutron stars. iRuffert fc Jankal (|l998l ) found that 
these shocks seriously pollute the direct collision vicinity 
with large amounts of baryonic material, which led them to 
rule out direct neutron star collisions as GRB progenitors. 
However, prospects for producing a GRB may be enhanced 
for off-centre collisions, which largely avoid these very strong 



suggest continued energy injection. A promising mechanism 
to power the observed la te-time activity is the fallback ac - 
cretion of the NS debris (|Rosswodl2007tlFaber et al.ll2006l ). 
and this effect could be en hanced in coUisional encounters 
(|Lee fc Ramirez-Ruiz|[2007l ). 

Collisions and close encounters between compact ob- 
jects occur predominantly in regions of h igh stellar density , 
such as the cores of globular clusters (e.g. iHills fc Davlll976l : 
lLeonard|[l989l : iFreitag fc Beii^bOOSh , or the central clusters 
of massive compact remnants that are expected in the cen- 
tres of galaxies (e.g. lAlexanderl2005l ). This pr oble m has been 
addressed in a number of pr evious studies (iKat z fc Canell 
ll996l : lHansen fc Muralilll998l ). In this paper we investigate 
an alternative possibility of compact object collisions in- 
duced by supernovae in binary star systems. We consider 
binary systems initially composed of two massive (>8M0) 
main sequence (MS) stars. The primary, i. e. the more mas- 
sive star, leaves the MS and explodes as a supernova (SN), 
giving birth to the first compact object, either a neutron 
star or a black hole. If the binary system does not get dis- 
rupted, we are left with a massive MS star in orbit around 
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a NS (or a BH). Since the primary loses most of its enve- 
lope mass prior to the supernova and the secondary is itself 
a massive star, there is a h igh probability for the binar; 



to survive the explosion (e.g. Ivaii den Heuvel fc Heisi 



! pmary 
s3 ll972l : 



IPortegies Zwart &: VerbuntlllQOd ). The life of the secondary 
star will soon come to an end. Low mass He stars in wide or- 
bits (Por6>0.25d) and those more massive th an ^^4 Mq do 
not so through common envelope evolution (jPewi fc Pold 
[2003) and the system will undergo the second SN explosion. 
At this point either a) the binary disrupts, or b) survives as 
a NS-NS or a NS-BH binary, or c) the newly formed NS is 
placed in a nearly parabolic orbit with a pericenter distance 
comparable to its tidal radius. In the last case we expect the 
NS to tidally interact with the companion compact object 
and rapidly spiral into it. This is the scenario considered 
here. 

Following the second S N, a debris disc c ould form 
around the newborn NS jZhang et al.l I2OO8I : iLin et al] 
I1991I : IChevaiieilll989l ). Such discs are thought t o be short 
lived systems (Trfj.^<.10^ vrs: IWang et aL l l2006bl : lEk si et al.l 
I2OO5I ). therefore they have no influence on NS-NS mergers 
(Tmerge~Gyr) Or coUisious in globular clusters. However, in 
our case the two compact objects will undergo a close en- 
counter within days of the explosion and the presence of 
a fallback di sc may have a significant effect on the event 
(|Popovll200^ ). 

The paper is organized as follows: in §[5] we investigate 
whether collisions induced by an asymmetric SN are feasible 
and how they depend on the binary parameters (eccentricity, 
orbital period, pre-SN mass, SN kick); in §[3] we derive the 
delay time between the precursor SN and the NSs encounter. 
The presence of a debris disc around the newborn NS is 
addressed in § 2] We discuss our results and summarize our 
conclusions in § [5] 



2 PROBABILITY OF A COLLISION 

We consider a binary system which already experienced a 
first SN explosion, and make two assumptions: 1) the sys- 
tem survived the explosion and now consists of a compact 
remnant and the remaining massive star and 2) it is in a 
circular orbit prior to the second supernova. 

The geometry of the system is illustrated in Fig. [T] we 
adopted a reference frame centered on the exploding star of 
mass Ml, and where the NS, of mass Mms, is at rest. The 
star is moving in a circular orbit, with separation ao and 
velocity Vorb = (GMo/ao)^^^, where G is the gravitational 
constant, and Mo = M\-\-Mms the total mass of the system. 

The direction of the kick is specified by the two angles 6 
and <j). As shown in Fig. [1] we defined 6 as the angle between 
the pre-SN orbital velocity Vorb and the kick velocity Vk, 
imparted by the SN explosion. The angle (j) is defined as 
the angle between the pre-SN orbital plane and the plane 
containing Vorb and Vk, where = (<^ = tt) if Vk lies 
in the pre-SN orbital plane and points outward (toward) 
the NS, and (j) = 7r/2 {(j) — 37r/2) if the plane containing 
Vorb and Vk is perpendicular to the pre-SN orbital plane 
and Vk has a Cartesian component parallel (anti-parallel) 
to the pre-SN angular momentum. If the kick velocity has 
no preferential orientation, then the values of cos 6 and 
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Figure 1. Geometry of the system: the reference frame is centered 
on the exploding star of mass Mi. The NS is at rest, while the 
star is moving in a circular orbit, with separation ao ^^ud relative 
velocity Vorb- The pre-SN orbit lies in the x-y plane, perpendic- 
ular to the plane of the page. The angle 9 is defined as the angle 
between the pre-SN orbital velocity Voj-b and the kick velocity Vk, 
imparted by the SN explosion. The angle </> is defined as the angle 
between the pre-SN orbital plane and the plane containing Vorb 
and Vk- 



have an equal probability of lying in the range [-1,1] and 
[0,27r], respectively. 

To study the effects of the second SN, we made the stan- 
dard assumption that the SN mass loss and kick are instan- 
taneous (i.e. on a timescale shorter than the orbital period), 
and neglected the impact of the expanding SN shell on the 
compact object companion, which have little or no effect 
on th e binary evolution (Kalogcra 199^; iFrvxell fc ArnettI 
Il98ll ). The specific angular momentum of the system is given 
by: 



|r X VI 



[GMf a(l 



,1/2 



(1) 



where r=(0, ao, 0) soon after the SN and V — Vorb + Vk 
is the resultant velocity of the newborn compact object. On 
the right side, Alf is the final mass of the system, a and e the 
post-SN orbital separation and eccentricity, respectively. In 
the following we consider the simplest case of two identical 
neutron stars (same mass and radius). 

In order to induce a collision the post-explosion or- 
bit must have a pericenter radius Rp=a{l — e) com- 
parable to the neutron star radius _Rjvs«10 km << ao. 
Therefore, only nearly parabolic orbits with eccentricity 
e > 1 — 2Rns /ao ~ 1 niay lead to such an encounter. From 
Eq.[T]it follows that the kick velocity must be such to signif- 
icantly reduce the orbital angular momentum of the system. 
This may happen a) if the kick velocity is much higher than 
the pre-SN orbital velocity, so that V ~ Vk, and directed in 
the radial direction; or b) if the kick velocity is comparable 
in magnitude with the pre-SN orbital velocity but directed 
in the opposite direction, so that V 0. 

In the former case (14 >> Vorb) the kick will unbind the 
binary, and the newborn NS will collide with the companion 
only if the kick direction points exactly toward it. The prob- 
ability of such alignment is small, Pcoii ~ Rns/o-o <10~**. 
Furthermore, such high values of Vk are very unlikely to oc- 
cur as they would imply an unfeasibly large population of 
isolated neutron stars for every observed binary pulsar. 

A collision between the two compact objects takes place 
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only if the following constraints on 9 and derived from 
Eq. [1] are fulfilled: 



( 



(2) 



}^(l + e)^cose^-^(l-e) (3) 



where 



Mo ao J 



(4) 



It follows immediately from Eq. |3] that if Vk < Vorb 
a close encounter is not possible, since cos 9 —1, and 
that the minimum kick magnitude to induce a collision is 
Vk = Vorb{i — ~ Vorb- FoT isotropic kicks, we derived the 
probability of a collision, pcoii, as a function of Vk by numer- 
ically integrating the following expression: 



Pcoll(Vfe) = - 



<i>MAx(Vk,0) 



d<f> (5) 



de sin e 

'OMAx(Vk) 

where the integration limits on <jf> and 6 are set by Eq. [5] 
and m The result is shown in Fig. [J] The values of Pcoii, 
quoted on the right y-axis, have been calculated for a kick 
magnitude Vfc=200 kms~^, a pre-SN mass Mi=4 Moand a 
NS mass Mns=1AMq. 

The probability shows a narrow peak at Vk/Vorb=i with 
a maximum value of Pcoii^S^'^/'^, then decreases rapidly. At 
Vk/Vorb=UMAx, with umax={1 + 2M} /MoY^'^ , the prob- 
ability drops by a factor of two. In fact, for a given mass 
loss, highly eccentric post-SN orbits are still bound only if 
Vk/Vorb < Umax, while higher values of Vk disrupt the sys- 
tem. In the latter case the two compact objects collide only 
if Vk has a Cartesian component directed toward the com- 
panion NS, i. e. only if 7r/2 < (j) < 37r/2. 

If the first born compact object is a BH, rather than 
a NS, the pericenter radius Rp must be comparable to the 
NS tidal radius, Rud^RNs{MBH /MnsY'^ , in order for the 
two objects to interact. In addition, the condition Rp>Rg, 
where Rg is the BH gravitational radius, avoids the possi- 
bility of the NS being swallowed directly by the BH. The 
latter constraint excludes BHs more massive than ~1OM0, 
while for less massive BHs the probability would change by 
a factor « {Mbh/Mns^^'^ - {Mbh/5 M^f'^ for a similar 
value of mass loss. 



2.1 Natal kicks 

Here we assume that kick velocities follow the birth veloc- 
ity distribution of isolated radio pulsars. The probability of 
NS-NS collision is therefore given by the convolution of the 
probability Pcoii(Vfe), derived in Eq.[Sl with the kick velocity 
distribution f{Vk): 



f 

Jo 



Pcm{V -Vk)f{Vk)dVk 



(6) 



lArzoumanian et all tOOt i modeled the kick distri- 
bution as two Gaussian components having characteris- 
tic dispersions (T;~90kms~^ and (t;i~500 km s~^. From a 




Figure 2. Probability of a collision as a function of the ratio 
Vk/Vorb between the kick velocity and the pre-SN orbital velocity. 
For a given mass loss, the final orbit is bound if Vk /Vorb < ^MAX- 
The values of Pcoii reported on the right y-axis have been cal- 
culated for a kick magnitude Vj,=200 kms^^, a prc-SN mass 
Mi=4M0and a NS mass of 1.4M0. 



larger sample of 73 young (<3Myr) pulsars, iHobbs et al.l 
(.2005) derived a mean three dimensional birth velocity of 
~400kms~'^ with a characteristic dispersion (j'~265 kms~^. 
The expression in Eq. |S] has been nurnerically integrated 
both for a bimodal distribution, as in Arzoumanian ct al] 
||2002|) and a broad Maxwellian distribution, as in 
iHobbs et all l|2005l ). Results are shown in Fig. [31 which re- 
ports the probability of a NS-NS collision as a function of 
the pre-SN orbital period. Calculations have been performed 
for a pre-SN mass of 4M0 (solid lines) and 10 M0 (dashed 
lines) . 

In order to obtain the distribution of post-SN orbital 
characteristics (velocity, separation, eccentricity), we ran 
Monte-Carlo simulations of a SN explosion in a circular 
binary system. We derived the effects of a SN on the dy- 
namics of a binary acco rding to the analytic formulation of 
iTauris fc Takend ( 19981 ) and the kick direction and magni- 
tude by Monte-Carlo techniques. We assumed an isotropic 
distribution of kicks {P^ — l/2n, Pg=sm9 /2) and draw the 
kick magnitude both from the bimodal and the Maxwell 
distribution of speeds used in our previous calculations. For 
each binary system, we ran 10* simulations and estimated 
the probability of a collision as the fraction of simulated 
SN explosions which lead to an encounter [Rp < 2Rns)- 
In Fig. [3] the results of the simulations are shown as filled 
squares (Maxwell distribution) and circles (bimodal distri- 
bution) overlaid on the solution of Eq. [6] The good agree- 
ment between the two confirms our analytic derivation. 

Throughout this work we assumed that kicks are 
isotropically distributed and uncorrelated with the binary 
properties, however this is a major source of uncertainty. 
Several physical mechanisms have been proposed to explain 
observ ations of high velocity pulsars (jLai 2001; Wang et al] 
l2006al . and references therein) and, for instance, some of 
them invoke a coupling between the dir ection of the kick a nd 
the stellar spin prior to the explosion (|Wang et al.ll2007l ). 

Fig. [3] shows that the probability of a collision event is 
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1 10 100 

Porb [d] 

Figure 3. Probability of a NS-NS collision as a function of the 
pre-SN orbital period Porb^ calculated for a pre-SN He star mass 
of 4M0 (solid lines) and 10 (dashed lines). The NS mass 
is 1.4 M0. Symbols report the results of numerical simulations: 
we assumed an isotropic distribution of kick directions, and ran 
simulations both for a bimodal distribution (red circles) and a 
MaxwcUian distribution (blue squares) of kick speeds. We re- 
stricted ourselves to the case in which the binary orbit is circular 
prior to the second SN. 



bound, it will cover almost the entire orbit before the en- 
counter; in this case the delay will be ~ 27ry^^^j — Af. 

For hyperbolic orbits (e > 1), we consider only those 
cases in which the NS is launched toward the companion 
((^~7r). The elapsed time since the explosion is: 

^'^yp = \/OTl"{[l + ^ + ^(^)' + 2^-(-"-l)]} (8) 

Fig. |4] shows the cumulative distribution of delays (left 
panel) and penetration factors I3=2Rn s / Rp (right panel) for 
the simulations described in § 12.11 The pre-SN mass of the 
He star was fixed to 4 M© , but the distributions presented in 
Fig.|4]do not depend on this particular choice. In fact the SN 
mass loss mainly affects the normalization of the probability 
curve rather than its shape (see Fig. |3]). 

As shown in Fig. U (left panel) , most of the collisions 
happen a few days after the SN. Though the events are 
nearly simultaneous, the SN will be detectable only a couple 
of weeks later, when it reaches its maximum optical luminos- 
ity (shadowed area). Only a small fraction (<5%) of events 
happen at later times (At >30 d). Off-axis collisions, which 
could more easily avoid a large baryonic contamination, are 
dominant (~80% with (3 < 2; Fig. g] right panel), while 
head-on collisions {13 » 1) are rare events (<5%). 



higher in close binary systems {Porb<W d). In such close sys- 
tems the presence of a binary companion can profoundly al- 
ter stellar evolution. It seems conceivable that it may some- 
how affect the SN ex plosion a nd the kick imparted to the 
stella r remnant (e.g. iPfahl et al.li200^ : IPodsiadlowski et al.l 
I2OO4I ). Tidal distortion, aspherical stellar winds, fast rota- 
tion, which in close binaries tends to be aligned to the or- 
bital spin, cause asymmetries in the star density and velocity 
distributions. The coupling between the stellar envelope and 
its core during the explosion is not understood yet, but one 
might expect that it can lead to an asymmetric mass ejec- 
tion, hence imparting a recoil velocity which is correlated 
with the orbital motion. 

The probability of an encounter would significantly 
change according to the kick geometry. If natal kicks are 
preferentially in the direction of the pre-SN orbital angular 
momentum, then the close encounters here described would 
not happen. By contrast, if the preferential direction is along 
the orbital plane in a cone with semi-aperture Aa, the num- 
ber of collisions would increase by a factor ~ 7r/2Aa< ^^^^ . 



3 DELAY SINCE THE SUPERNOVA 

These collisions are triggered by an asymmetric SN explo- 
sion and will happen soon after it. The time interval At that 
elapses since the SN can be derived from the known orbital 
parameters: 



(7) 



where the post-explosion eccentricity e and separation 
a have b een calculated for each s imulated SN following the 
recipes of lTauris fc Takend l|l998l ). If the newborn NS is shot 
in the direction opposite to the companion (<^~0) but still 



4 THE EFFECT OF A FALLBACK DISC 

In a SN some of the mass is expected to remain bound to the 
central object and carry suff i cient angular momentum to set - 
tle into a disc (iMichellllQSSi : iLin et al.lll99ll : IWooslevlll993l ). 
The presence of such discs around NSs has been invoked 
to explain the peculiar class of anomalous X-ray pulsars 
(AXPs) and indeed the first observational evidence of their 
existence has been d iscovered around the AXP 4U 0142-1-61 
dWang et al.ll2006bl). Yet it is unclear whether this disc is 
passive (IWang et al]l2006d) or affect s the propert es of the 
newly formed NS l|Ertan et al.ll2007l ). Instead, if a fallback 
disc persists after a SN in a binary system, as considered 
here, it may drastically increase the cross-section of an in- 
teraction event. 

Numerical simulations of stellar evolution find that 
the specific angular momentum distribution ranges from 
j ~10^''-10^'' cm^ s"^ (in the stellar core) to j >10^^^cniV_i 
(in t he outer layers) prior to the explosion l|Heger et al.l 
l2005h . The supernova shocks traversing the star likely in- 
duce an extensive mixing between the stellar layers and 
an appreciable fraction of fallback material may have an 
angular momentum as high as j '^lO^'^ cm^ s~^. By as- 
suming that conservation of angular momentum holds dur- 
ing the collapse, this sets the initial size of the fallback 
disc Rns < Rd < 10^ cm. For typical v alues of natal kicks 
(<1000 kms~^ [Chatteriee et al1l2005l ). the disc is not left 
behind, but remains bound to the new born NS. The system 
is therefore composed of a NS (or a BH, formed in the first 
SN explosion) and a new born NS (formed in the second SN) 
surrounded by a fallback disc. 

After an initial transient phase the disc spreads to larger 
radii as -Rd oc (t /ri,)^^^, where Ti, »10^-10'*s is the local vis- 
cous t imescale (|Lvnden-Bell fc Pringldl 19741 : ICannizzo et al.l 
Il990l :l Menou et al.ll200ll ). and reaches a size of ~10^ cm at 
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Figure 4. Left panel: cumulative distribution of delay times since the SN explosion; the shadowed region indicate the time interval of 
the SN peak . Right panel: cumulative distribution of penetration factors, /3, defined as the ratio between the sum of the NSs radii and 
the post-SN pericenter distance, /3=2R]vs/-Rp. 



30 d since the SN. The chance for the disc to interact with 
the stellar companion shortly after the second SN is there- 
fore cx {Rd/ Rns)^^'^ 100 times higher than the probability 
of a direct collision derived in §(2] The perturbation induced 
by the passage of the first NS through the disc extends up 
to the NS accretion radius = GMns /v^ei ~ Rd, where 
Vrei is the relative velocity between the disc and the first 
NS. This is comparable with the disc radius itself and there- 
fore most of the mass in the d i sc (10'^^ Mp , < Md < . 1 M© ; 
IChevaheil [l989l : | Lin et ahl Il99ll : iFrverl l2009l i is transferred 
to the NS on a dynamical timescale of ~2 {Rd/10^ cm)^'^^ s. 
This episode of accretion may result in a luminous X-ray 
outburst a few days/weeks after a SN. The precise nature of 
this outburst depends on the details of the interaction and 
the associated gas dynamics but it is highly likely to involve 
accretion rates close to Eddington. 



The presence of a debris disc in a binary system 
may also lead in some cases to a rapid shrinkage of 



5 SUMMARY AND DISCUSSION 

We showed in § 2 that a supernova may be followed by the 
collision of two compact objects. A close encounter requires 
the explosion to be asymmetric (spherical explosion always 
led to an expansion of the orbit) , and happens only for a nar- 
row range of kick magnitudes and directions, namely when 
the kick velocity is comparable to the pre-SN orbital ve- 
locity and opposite in direction. These events are therefore 
expected to be rare. 

Similarly to mergers, collisions of compact objects 
end up forming an accr etion disc/black hole system 
(|Lee fc Ramirez-Ruid lioOTl ). as commonly invoked as the 
central engine powering GRBs. One of the most important 
open issue related to (all) the GRB central eng ines is the 
"baryonic contamination" problem (|Piranlll999r and refer- 
ences therein): How can a violent explosive event involv- 
ing a few solar masses of baryons channel a huge amount 
of energy into a region that is almost completely baryon- 
depleted? The popular answer to this question in the con- 
text of neutron star mergers has been the centrifugally evac- 
uated funnel regio n that forms above the poles of the cen- 



the o rbit, and speed up 



the merger of the two ob- tral remnant (e.g . Davies et al.l 1 19941 : iRosswog et al. I I2OOOI : 



iects (iGoldreich fc Tremaindll980l : iLin fc Papaloizoulll986l : 
lArmitage fc NataraianI |200^ ) . For fallback discs of mass 
Md<<O.lM0, the characteristic time scale of this proces s 
is Tshrir^k/r^, > {Md + MNs)/Md » 1 ijLodato et a l.l l2009h . 
The effects of such an interaction are not likely to be impor- 
tant since, as pointed out above, most of most of the disc 
mass is transferred to the NS companion on a much shorter 
timescale. 



iRuffert fc Jankall200ll ). ahhough it had been realized rela- 
tively early on that this loophole could possibly be threat- 
ened by the ablation of baryonic m aterial via th e extreme 
neutrino luminosities of ~ 10^'^ e rg/s (|Ruffert fc Ja nka 199^; 
IRosswog fc Ramirez- Ruij|2002l ). The first calculations that 
accou nt for neutrino heat ing processes in a NS merger rem- 
nant (|Dessart et al.l [20091 ) seriously downsize this loophole. 
For at least as long as the central, super-massive NS has not 
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collapsed into a black hole, the remnant drives a very strong 
baryonic wind (Mu,~ 10~^ Mqs~^) into exactly the above 
mentioned funnel region, making the launch of a relativistic 
outflow impossible. This conclusion may need modification 
once a black hole forms. 

For the case of neutron star collisions the situation is 
even less promising. At least head-on collisions (/3>>1) with 
their violent shocks heavily spoil their surroundings with 
baryonic material even if neut rino-driven winds are disre- 
garded (|Ruffert fc JankallT997l ). However, we have shown 
that only a small fraction (<5%) of collisions are expected 
to be near-central (cf. Fig. IH right panel), while most of 
them off-centre. In this case the outcome may be more fa- 
vorable to GRBs. Off-centre collisions lead to much weaker 
shocks, and a much smoother subsequent coalescence, likely 
dispersing a substantial amount of debris material in highly 
eccentric orbits. This material in turn can pr oduce late-time 
X-ray flares simialr to tho se seen in GRBs (|Rosswogll2007l : 
iLee Ramirez-Ruij2007l '). Collisions between a stellar mass 
BIf and a NS are even more likely to produce a short GRB 
since they avoid two problems from the beginning: there will 
only be moderate shocks in the disk formation phase after 
the disruption and there is never a hot, massive NS as a 
major source for the neutrino wind. So in the latter cases 
it may be possible to drive the ultra-relativistic outflow to 
power a short-duration GRB. 

In the case of NS collisions induced by a SN explosion, 
for a GRB to occur the emerging jet must have a mini- 
mum power to break out of the remnant layers maint ain- 
ing a high Lorentz facto r (jMacFadvcn ot al. 2001; Mat zneJ 
l2003l : IZhang et al.|[2003 '). This condition is particularly re- 
strictive when the delay time between the SN and the col- 
lision event is short (t<<l d). As the supernova remnant 
(SNR) expands, the requirements on the jet energy and 
beaming angle decrease, therefore even low luminosity and 
less collimated jets could pass this constraint. Furthermore, 
if the SNR radius is larger than the radius for internal shocks 
{Ris ~10^^-10^^ cm), the prompt GRB emission can plau- 
sibly occur within it, as the precursor SN has the advantage 
to clean the circumbu rst environment from baryon pollution 
(|Vietri fc Stellalll998D . thus favoring a relativistic expansion 
of the GRB jet within the remnant. A major problem for 
longer time delays is whether the GRB emission can be ob- 
served or it is obscured by the dense SNR shell. In fact, for a 
radial velocity of the ejecta Vej~10^ kms~^, the SNR radius 
is i? ~ 10^'* cm after 1 d, and the resulting column density 
is nii~102**(M,heH/Mo)cm-^ 

Both these problems (the jet break-out and the SNR 
optical thickness) are not present if the geometry of the 
SNR shell is not spherical , but anisotropic either because of 
an as ymmetric explosion ([Leonard et al.|[2006l : iMaeda et al.l 
[2008') or because subject to hydrodynamic instabilities 
(Chevalier ot al. 1992; Blondin & Ellison 2001). The SNR 
might not interfer with the jet and the GRB emission may 
escape from underdense regions. In this case, a short GRB 
(or an X-ray outburst, see §|4]) could follow the SN explo- 
sion. In contrast to long G RBs, connected to the peculiar 
broad line Type Ic SNe (jWooslev fc Bloom i2006 ') , bursts 
originated through this channel will be associated with any 
type of core-collapse SN. The supernova always preceeds the 
NSs close encounter, but according to the distribution of de- 
lay times shown in Fig. [4] the high energy transient should 



be observed first, while the SN emission will be visible a few 
days later. For small kick velocities a longer delay between 
the SN and the NSs collision is also possible: ~10% of our 
events happen a couple of weeks after the stellar collapse, 
and a small fraction, <5%, happens more than one month 
later. 

According to the proposed scenario, double NSs and 
close encounters have the same progenitor systems, i.e. a 
massive binary system which undergoes two SN explosions. 
A fraction of NS-NS binaries will mer ge within a Hubbl e 
timescale, giving birth to a short GRB (|Eichler et al.|[l989l ). 
By assuming that the main progenitors of short GRBs are 
binary mergers (either NS-NS or NS-BH), the observed rate 
of collisions n°oa is simply proportional to the rate of short 
GRBs hsGRBs- 

• obs r ■ ^coll r . 

n-cou = Je^ricou = JOm nsGRBs (9) 

nmergers 

where /g^ and /e„ are the beaming corrections for bursts 
produced through collisions and mergers, respectively; 
ncoii/nmergers is the relative strength of the two channels, 
the former giving rise to prompt collisions and the latter 
le ading to merging NS s . Foll owing the method described 
in iKalogera fc Lorimej (|200(]| ) we derived an upper limit 
of ncoii/nmergers'^W~^ ■ With a fuU-sky rate of ~170 short 
GRBs per year, as derived from the BATSE GRB Cata- 
logu^il, and by assuming the same beaming factor of short 
GRBs originated from double NSs {fg^ =/e„ ) , we estimated 
for our events an observed rate of <0.2 yr~^. This num- 
ber increases by more than a factor of 10'' when considering 
the interaction with a disc of supernova debris, which may 
produce a luminous X-ray outburst. 

A wealth of new transient phenomena are to be dis- 
covered in the near future, when the sky will be con- 
tinously monitored both at high energies (e.g. Swtft, 
Fermi, MAXI) and at optical/NIR wavelenghts (e.g. PAN- 
STARRS, LSST). CoUisions induced by a SN explosion, 
though rare, might be detected by forthcoming surveys. Such 
events will be located in young stellar environments, such 
as stellar star clusters or star forming galaxies, or likely 
at redshift z > 1, when the star formation rate had its 
maximum (e.g. lMadau fc Pozzettl2000l : lHopk ins fc BeacomI 
12006! ). Though so far short GRBs have been found at 
moderately low redshifts < z >^ 0.5 and in regions of 
low/moderate star formation, we note that in the growing 
sample of GRBs with known redshifts a few bursts, charac- 
terized b^;_aj;atlMr_shor^jn duration (e.g. 090426 at 
z=2.6: ILevesgue et al.ll2"009^ ■ may have originated from the 
direct collisions described in this paper. 
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